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Abstract. We present a theoretical study of bose condensation and specific heat of non-interacting bosons 
in finite lattices in harmonic potentials in one, two, and three dimensions. We numerically diagonalize 
the Hamiltonian to obtain the energy levels of the systems. Using the energy levels thus obtained, we 
investigate the temperature dependence, dimensionality effects, lattice size dependence, and evolution to 
the bulk limit of the condensate fraction and the specific heat. Some preliminary results on the specific 
heat of fermions in optical lattices are also presented. The results obtained are contextualized within the 
current experimental and theoretical scenario. 

PACS. 3.75.Lm, 03.75.Nt, 03.75.Hh 



I 

o . 

^ , tical lattices thus far mainly concentrated on quantum 



1 Introduction [SlI^ITllSllWnTOfT^IT^IWT^in^lfTTllTH] of bosons in op- 



^ , Bosons in optical lattices provides a microscopic labora- 



phase transitions between various possible phases such 



tory for exploration of the properties of many-boson sys- i, ^ • , , , • , , , 

$H as Mott msulator phase, density-wave phase, and bose- 

d ' 

■ - - ' tems with unprecedented control on number of bosons per i i i ■ ■ r ^ 

condensed phase, h mite temperature properties of bosons 

in the combined optical lattice potential and the confin- 
ing harmonic potential have not received much attention 
except in some studies [T91I20] where, again, the focus is 
on Mott transition. 



site, boson kinetic energy, and boson-boson interaction 
strength [T1[51[31I311S] . Experimentalists have explored [Tl[2j 
[3ll4] quantum phase transitions, excitation spectra, con- 
densate fraction as a function of boson-boson interaction 
strength, and properties of boson Tonks gas in one, two, 
and three dimensional optical lattices. Theoretical studies 
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In this paper, we present a theoretical study of bosons 
in optical lattices in a harmonic potential in one, two, and 
three dimensional optical lattices. In particular, we study 
the temperature dependence of the condensate fraction 
and specific heat. This study is useful for the following 
reasons. After the discovery [^TU^^ll^ of Bose-Einstein 
condensation in harmonic traps (for reviews see Refs. [241 
[25 l [26 l l27] ) ■ one of the important studies done was to mea- 
sure the temperature dependence of the condensate frac- 
tion and its comparison with theory, and such measure- 
ments for bosons in optical lattices can be expected in 
near future. Thermodynamic properties of bosons in opti- 
cal lattices are of considerable interest in their own right. 
Though no measurement of specific heat of bosons in a 
harmonic potential and of bosons in a lattice with har- 
monic potential have been reported thus far, such mea- 
surements have already been reported [28] for fermions in 
harmonic potential, and it is not unreasonable to expect 
such measurements for bosons as well. Furthermore, ex- 
periments are performed on bosons in finite lattices, and so 
it would be useful to understand the properties of bosons 
in finite lattices to see finite size effects and their evolution 
towards the bulk limit. The motivations and relevance of 
the studies presented in this paper are, hopefully, clear to 
the reader from the foregoing. 

As stated, we present studies of lattice bosons in har- 
monic traps. We have considered non-interacting bosons 
in a periodic lattice in one-dimension (Id), a square lat- 
tice in two-dimensions (2d), and a cubic lattice in three- 
dimensions (3d) . The bosons are considered to be moving 
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in the combined lattice and harmonic potentials. We nu- 
merically diagonalizc the Hamiltonian matrix for bosons 
in finite lattices in a harmonic potential, and concentrate 
on calculations of the temperature dependence of ground 
state occupancy, specific- heat, finite-size effects, and di- 
mensionality effects. This work is presented in the Sec- 
tions 2 and 3. In Section 4, some preliminary results on 
our recent work on fermions in optical lattices in a har- 
monic potential is presented, and conclusions are given in 
Section 5. 

2 Model and method 

In this Section, we will describe the model and the method 
followed. Consider non-interacting bosons under the com- 
bined influence of a periodic optical lattice potential and 
an overall harmonic confining potential. The Hamiltonian 
of this system is, 

<ij> i i 

where t is the kinetic energy gain when a boson hop from 
site i to its nearest neighbor site j in the optical lattice, 
c\ is the boson creation operator, K (> 0) the strength 
of the harmonic confining potential, the locator of the 
site i with respect to the origin which is at the center of 
the lattice, Ui — c\ci the boson number operator, and ^ 
the chemical potential. For bosons under the influence of 
only the harmonic potential bose condensation is by now 
well understood [Ml[^51l^51l27j . Our aim is to consider the 
case of finite t and finite K . After writing the Hamiltonian 
(Eq.l) in the single particle site basis, we numerically di- 
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agonalize it for various values of K and for various lattice chemical potential and the population in the lowest level 

sizes in Id, 2d, and 3d. The energy levels thus obtained (A^o) are determined using Eqs. (2) and (3). 

are used in calculations of the ground state occupancy and The total energy of the system of bosons is given by 
specific heat as a function of temperature. We will see that 

Etot = Y,N{Ei)E,. (4) 

any non-zero value of K leads to substantial changes in i=o 

both these properties. We have calculated the energy lev- The temperature derivative of Etot gives the specific heat 



els and the bosonic properties using open boundary con- 
ditions. It may be noted here that for the Id case, the 
eigenfunctions are eigenenergies may be exactly obtained 



hh 



lo 



(5) 



where, 



in terms of the Mathicu functions as shown by Rey and j _ p(^E- — ^) (6) 

j— 

collaborators |29j . They have given approximate analyt- ^ 

h^Y.E,F{E,-ti), (7) 

/2 = (8) 

z^O 

/3 = (9) 



ical expressions for the eigenenergies and eigenfunctions. 
We compared our eigenenergies for Id systems with their 
approximate result (Eq. 15 in their paper) in the tunnel- 
ing dominated regime and found that their eigenenergies 
are very close to our exact results for large values oit/k. 

The ground state occupancy and the specific heat are i3(E --ti) 

F{E, - ^) = ^ 2 . (10) 

calculated in the following way. The population N{Ei) of ^el^(Ei-fj.) _ 



z=0 

In the above equations 



a state with energy Ei is given by the Bose distribution 



^(^^) - e^iE^^-l ' (2) 



3 Results and Discussions on bosons 



Before presenting the results of our calculations, some 

where /3 = l/ksT and ks is the Boltzmann constant, general remarks are in order. In the absence of a har- 

The chemical potential is determined from the number monic confining potential, the bosons in a periodic lattice 

equation undergo a bose condensation phase transition only in 3d 

Jzz. where the boson Density Of States (DOS) in the thermo- 
iV = ^7V(i?,), (3) 

»=o dynamic limit vanishes at the bottom of the band. The 

where Ei with i— 0, 1, 2 , im denotes the energy levels phase transition is not possible for bosons in Id or 2d pe- 

from the lowest to the highest level of a boson in a finite riodic lattices since the DOS at the bottom of the band 

lattice in the presence of the harmonic confining potential, either diverges or remains finite, and the condition Nq — 

For a given number of bosons, at any temperature, the and ^ = Eq is not satisfied at any temperature. Here, Eq 
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is the energy of the lowest level or bottom of a band. For 
free bosons in a harmonic confining potential, the phase 
transition is possible in 2d but not in Id [3D]. For finite 
size systems and for finite number of particles there is no 
true phase transition using the above criteria. However, 
at sufficiently low temperatures, a macroscopic number of 
particles will occupy the ground state and thus bose con- 
densation will occur [31j . With these general remarks, we 
go to our results. 

3.1 One Dimensional case 

In Fig. 1, the variation of the condensate fraction (Nq/N) 
with temperature is shown for 1000 bosons in a Id lat- 
tice (of lattice constant a) for k — Ka^ — 0.01 (in the 
energy scale of t 1) for different lattice sizes (N/). The 
occupancy in the lowest level is considerable even at high 
temperatures for smaller size systems and it decreases with 
increasing size. For k = 0.01, no size effect is seen for N; > 
400. We have presented a curve for N; — 1000 which may 
be considered as an infinite size result for this value of k. 

In Fig 2. we have shown the temperature dependence 
of the condensate fraction (Nq/N) for a one dimensional 
lattice of size TV; = 1000 in presence of harmonic potential 
of strength k — 0.01 for different numbers of bosons. The 
Nq/N is plotted against T/Tq. Here, Tq is determined by 
setting A'o = and fi — Eq in the number equation (Eq. 
3) to obtain 
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Note that in the thermodynamic limit and when the DOS 
vanishes at the bottom of the boson energy band, Tq gives 
the phase transition temperature. For finite size lattices 
To may be considered as the condensation temperature. In 
the absence of a harmonic trap, Nq/N is appreciable for 
bosons in a Id lattice even at T = 2Tq (dash-dot curve 
in Fig. 2). It decreases rapidly with the application of 
the harmonic potential and also decreases with increasing 
number of bosons for a fixed value of k. For a bose gas in 
a harmonic trap (dashed curve in Fig. 2), Nq/N almost 
vanishes at high temperatures. In Fig. 3, the variation of 
To with number of bosons (N) is shown for different values 
of the harmonic potential. The To increases almost linearly 
with N. For a fixed N, it increases with increasing strength 
of the harmonic potential. 

In Fig. 4, we have plotted the specific heat per particle 
against T/Tq for different values of k for Ni = 1000 and 
N — 600. It is seen that the specific heat decreases with 
increasing value of k. For k — 0.001, a broad peak is seen 
in the specific heat, which is a signature of the finite size 
of the lattice. For such a small value of k, even a size of 
1000 sites is not enough to confine the bosons within the 
harmonic trap (i. e., some of the bosons may occupy the 
boundary sites), hence finite size effects would show up. 
As the strength of the harmonic potential increases, the 
number bosons at high-energy sites at the edges of the lat- 
tice is insignificant and the finite size effect disappears. We 
have also shown Cy for bosons under the influence of only 
the harmonic potential. One notices at once that while 
the dependence of C„ on the scaled temperature (T/To) 
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for free bosons in a harmonic potential is independent of 
k, the same is not true for lattice bosons. It is found that, 
switching on the lattice potential substantially reduces the 
Cy. We also note that reducing the strength of harmonic 
potential for the lattice bosons drives the system towards 
the limit of free bosons in a harmonic potential, which is 
counter intuitive to what one might naively expect. In the 
very low temperature range, all the curves are very close 
to each other. This can be understood if we look at the 
single particle spectrum of a large lattice. In this case, the 
low energy part of the band is approximately parabolic, 
and so one can transform the problem to a free particle 
in a harmonic potential. Hence one expects the proper- 
ties of lattice bosons in a harmonic trap will be close to 
that of free bosons in harmonic trap at sufficiently low 
temperatures. 

3.2 Two Dimensional case 

In this section we present our results on bosons in two 
dimensional optical lattices in a 2d harmonic potential. 
These results arc presented in Figs. 5-10. Fig. 5 shows the 
variation of the condensate fraction with temperature for 
different lattice sizes for A; = and 0.01, respectively. A 
macroscopic condensation occurs at a lower temperature 
for larger size of the system. The dashed lines in Fig. 5 
exhibit the effect of harmonic potential on lattice bosons 
in 2d. In the presence of a harmonic confining potential, 
No/N versus T curves become almost size independent 
for lattice sizes of 80 x 80 and larger, unlike the case of 
A; = (solid lines) where the condensation temperature 
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decreases rapidly with the size of the system. For largest 
lattice size shown, in presence of the harmonic potential 
our results may correspond to the results in the thermo- 
dynamic limit. It may be mentioned that for larger num- 
ber of bosons a larger lattice size is needed to reach the 
infinite size limit. The Tq is clearly seen to be consider- 
ably increased with the introduction of the harmonic po- 
tential. We also note that the condensate fraction grows 
faster with decreasing temperature in 2d as compared to 
that in Id (see Figs. 1 and 2). In Fig. 6, we have shown 
the variation of Tq with boson number. The figure shows 
that, for a fixed lattice size and A'^, Tq increases with in- 
creasing strength of the harmonic confinement potential. 
It also shows that the bosons under the combined effect 
of lattice and harmonic potentials have a higher transi- 
tion temperature compared to bosons with either of the 
potentials absent. In comparison with Id, Tq is found to 
increase rapidly for small N and a monotonic increase for 
larger N, clearly showing a dimensionality effect. 

In Fig. 7, we have exhibited the effect of lattice size 
on the specific heat of bosons for a fixed N and k. For 
small lattice sizes, we find a flat region in Cy for tempera- 
tures below Tq. With increasing lattice size the flat region 
evolves into a peak which becomes sharper as the size of 
the lattice approaches the bulk limit. If one compares the 
C„ results for 80 x 80, 100 x 100, 150 x 150 lattices, one 
finds that in the low temperature range the value of Cy is 
same for all the three curves, while in the high tempera- 
ture range Cy is smaller for smaller size systems. The latter 
is a signature of the size effect. For higher temperatures 
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more and more bosons would occupy the boundary sites gularity for a 2d square lattice. Any finite harmonic po- 
of the lattice for a given value of k leading to appreciable tential leads to drastic changes in DOS by wiping out the 
finite size effects. Given the great control experimentalists van-Hove singularity and fiattening the DOS in an energy 
have demonstrated, these effects may be within the reach range which increases with increasing harmonic potential 
of experimental observability in near future. strength (k). This range also increases with increasing size 

T i TiirTiri of the lattice (not shown in the figure). The DOS at the 

In order to get an understandmg ot the enect ot har- ^ ' 

, 1 i ii ii -n ^ i r 1 bottom of the band vanishes as the size of the lattice ap- 

monic potential strength on the specihc heat of lattice 

, , 1 • TT o /-( r i- rrn/m proachcs thc infinite limit, and in the low energy range 

bosons, we have shown in i ig. 8, C„ as a iunction ot J /Jo oj o 

r . , r I rni • J2 1 ii i ii ^ it shows almost linear increase with energy before taking 

tor various values of k. ihis fagure shows that the Cv 

r j-rr . i i i thc coustant value. The DOS for bosons in a 2d lattice in 

curves tor dinerent k values intercept at a temperature 

T/To « 1.1 for the boson number {N = 3840) consid- P'^^^^"^^ ^ harmonic trap may thus be approximately 
ered. Below this temperature Ci, is lower for higher k, S^"^^^ by 

while above it increasing k increases C^,. This is in con- p{E) = for < E < Ei 

trast to the behavior of bosons in Id (Fig. 4) where in- 
creasing k decreases except at very low temperatures. 
The peak in Cy curve is seen to decrease with increasing 
harmonic confinement. The behavior of Cy shown in Fig 
8 results from an interplay of the effects of the harmonic 
potential and finite size of the lattice. The cross-over ef- 
fect for T/Tq > 1.1 occurs because in this temperature 
range bosons are delocalized over the entire lattice so that 
they sense the existence of boundaries. This is clear from 



= -/El for Ei<E <E2, (12) 

where all the energies are measured from the bottom of the 
energy spectrum. It is found that 7 is strongly k depen- 
dent and decreases with increasing k, while Ei is almost 
k independent. For free bosons in a 2d harmonic trap, the 
DOS has a linear dependence on E in the entire energy 
range (0 < E' < 00). For p{E) oc i?"^^, the specific heat 
per boson for T < Tq is given by[57] 



Fig. 9 in which C„ is shown for a larger lattice (150 X 150) — ^ ^ q;(q; -|- 1) ^^^^ XTq) ' 

and one notices that the crossing has disappeared in the , >/n v^oo/T/„\-ii tt 

^ where = i^n=\\^i'^ ) Kiemann zeta tunction. 

temperature range shown. c--(-ru -u ■ ^ o/-<u 

^ ° Since for tree bosons m a harmonic trap a = 2, Ct, shows 

In order to get an insight into the dimensionality ef- a (T/Tq)^ behavior for T < Tq. This feature is seen in 
fects, we calculated the one boson density of states for our calculated results (the dotted curve in Fig. 9). For 
a lattice boson in a harmonic potential. The normalized the lattice bosons in a harmonic trap in the low tempera- 
DOS for the 2d case is shown in Fig. 10. In the absence ture region, only the low energy linear part of of the DOS 
of the harmonic potential, the DOS has a van-Hove sin- is occupied. Hence shows behavior similar to that of 
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free bosons in a harmonic trap. When the temperature is 
appreciable so that the flat portion of the DOS is also oc- 
cupied, the effective value of a decreases. This results in a 
slower increase of Cy with T/Tq and lower value of Cy com- 
pared to the free bosons case. These features are promi- 
nent in Fig. 9 for higher values of k which have higher Tq 
values. 

3.3 Three Dimensional case 

In Fig. 11-13, we have presented our results for bosons in 
3d cubic lattices in a harmonic potential. The tempera- 
ture dependence of the condensate fraction (Fig. 11) and 
the boson number dependence of the condensation tem- 
perature (Fig. 12) shows variations similar to those found 
for 2d systems. We find that Tq increases with increasing 
harmonic potential strength as well as with the total num- 
ber of bosons in the lattice. In Fig. 13, we have shown the 
effect of harmonic potential strength on the specific heat 
of 3d lattice bosons. The peak value in Cy at T/Tq « 1 
is lower for higher values of k. For the 50 x 50 x 50 lat- 
tice, Cy is lower for lower values of k at high temperatures 
{T/Tq > 1.3). This results from the finite size of the lat- 
tice, as mentioned previously. In comparison with a 2d 
system, the harmonic potential strength has a lesser ef- 
fect in 3d (see Figs. 9 and 13). In Fig. 13, we have also 
plotted the Cy for free bosons in a harmonic trap. The 
Cy shows a (T/Tq)^ dependence for T/Tq < 1, which is a 
known result [^7]. The figure shows that the Cy for lat- 
tice bosons in a harmonic trap has the same temperature 
dependence in the low temperature range. We have cal- 
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culated the DOS for the bosons in the 3d lattice (of size 
50 X 50 X 50) confined in a harmonic trap and found that 
it has an E'^ dependence in the low energy range. This 
accounts for the (T/Tq)^ behavior of the Cy for the lattice 
bosons in a harmonic trap in the low temperature range. 

In the work presented on bosons, the boson-boson in- 
teraction is neglected. We believe that our results would 
approximately hold when interaction energy (U) is much 
smaller than the hopping energy (t). In an optical lattice, 
the interaction energy (U) depends on the ratio between 
the depth of the optical potential to the recoil energy. A 
weakly interacting regime may be obtained by adjusting 
the value of the potential depth to a low value [Ill32|. In 
the weak interaction regime, the interaction induced de- 
pletion effects may not be significant. 

4 Results on Fermions 

Due to recent interest in fermions in optical lat- 

tices, we have studied single (spin) component fermi gas in 
optical lattices with harmonic confinement in Id, 2d, and 
3d. Another interest in studying this fermionic system is 
that the bosons in the strongly interacting regime behave 
similar to fermions to avoid double occupancy and mini- 
mize inter-particle repulsion energy and energy spectrum 
of strongly interacting bosons in Id is very close to that 
of the non- interacting fermions |29J. 

Our results for the specific heat for different values 
of the trapping strength and for different dimensions are 
presented in Figs. 14-16. The specific heat curves show a 
linear variation at low temperatures, which is a charac- 
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teristic of a degenerate fermi gas. At high temperatures, 
the specific heat shows a more or less flat behavior. The 
specific heat increases with decreasing k values except at 
low temperatures. Similar results are also noted for bosons 
(see Figs. 4 and 9). With increasing k value, the specific 
heat for fermions approaches the corresponding value for 
the k = case. 

We have determined the gradients of the specific heat 
curves at low temperatures. For 3d, this value is 9.76 for 
k = 0.01 and 4.98 for k = 0. The former is close to the 
corresponding value (tt^) for free fermions in a harmonic 
trap while the latter is close to that ('7r^/2) of free fermions. 
For the 2-d case, this slope is 7.2 for k = 0.01 and 3.37 for 
k = 0, while for 1-d it is 3.95 for k = 0.0001 and 1.69 for 
fc = 0. Wc find that in all dimensions, the values of the 
gradient of the specific heat curves at low temperatures 
for small [k/zt) (where, z is the coordination number of a 
lattice), are close to those of free fermions in a harmonic 
trap. The reason behind it lies in the fact that the single 
particle density of states for small k values in an optical 
lattice shows the same energy dependence as that of the 
free fermions (bosons) in a harmonic trap at low energies, 
as discussed in previous sections. 

5 Conclusions 

In this paper, we have presented a study of non-interacting 
bosons under the influence of combined optical lattice and 
harmonic potentials in one, two, and three dimensions. 
The condensate fraction in Id shows a faster reduction 
with increasing temperature in the low temperature range 
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compared to that in 2d and 3d. The condensation tem- 
perature in Id increases linearly with increasing number 
of lattice bosons. In comparison, Tq in 2d and 3d shows 
a fast growth for small number of bosons and a mono- 
tonic increase thereafter. It is found that, for a given num- 
ber of bosons, the condensation temperature is higher for 
bosons in the combined harmonic and optical potentials 
compared to the cases in which either of the potentials is 
absent. 

The specific heat results show several interesting fea- 
tures. In Id, specific heat of lattice bosons in harmonic 
potential is found to show a very slow growth with tem- 
perature, except in the low temperature range where the 
growth is relatively faster. In 2d and 3d, C„ of lattice 
bosons in a harmonic trap has a peak at the condensation 
temperature. When the lattice size effects are important in 
C„, we find that the Cy per boson is lower for smaller lat- 
tice sizes except at low temperatures. In 2d, fiat regions 
are observed in C„ for a substantial temperature range 
below the condensation temperature for relatively smaller 
lattices. When the lattice size effects arc negligible, we find 
that the lattice bosons in a harmonic trap has a consider- 
ably reduced specific heat compared to that for free bosons 
in a harmonic trap for temperatures around and above Tq. 
In all dimensions, the low temperature is very close to 
that of free bosons in a harmonic trap. The specific heat of 
the lattice bosons is strongly dependent on the strength of 
the harmonic potential in contrast to that of free bosons 
in a harmonic potential. In all dimensions for large size 
systems when the lattice size effects are negligible, the Cy 
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is lower for higher values of k around and above Tq. How- 4. M. Kohl, H. Moritz, T. Stoferle, C. Schori, and T. Esslinger, 

ever, the C„ curves for different k values cross each other J- Low Temp. Phys. 138, 635 (2005). 

above Tq when finite size effects are present. Considering 5. D. Jaksch, C. Bmder, J. I. Cirac, C. W. Gardiner, and P. 

recent interest in fermions in optical lattices, we presented Zoller, Phys. Rev. Lett. 81, 3108 (1998). 

some preliminary results on the specific heat of fermions 6. D. Van Oosten, P. van der Straten, and H. T. C. Stoof, 

in optical lattices in a harmonic potential in Id, 2d, and P^iys- R-ev. A 63, 054601 (2001). 

3d. The specific heat is found to show a linear tempera- 7. W. Zwerger, J. Opt. B: Quantum Semiclass. Opt. 5, S9 

ture dependence at low temperatures and a flat behavior (2003). 
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pure lattice fermion results. The temperature dependence (2003). 
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harmonic trap is governed by a complex interplay of the ^^"^^^^ (2003). 
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Fig. 1. The variation of condensate fraction with tempera- 
ture for 1000 bosons in Id finite lattices in a harmonic po- 
tential with k = 0.01. The size of the lattices are: N( = 30 
(dots), 50 (dashes), 100 (dash-dot), and 1000 (solid line). ksT 
is expressed in energy units of t = 1. 
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Fig. 2. The variation of condensate fraction with temperature. 
Solid lines are for bosons in a Id finite lattice of size 1000 in a 
harmonic potential with k = 0.01 for boson numbers N = 100 
(top), 1000 (middle), 10000 (bottom). The dashed fine is for 
10000 free bosons in a harmonic trap with k = 0.01. The dash- 
dot line is for 10000 bosons in a Id lattice of 1000 sites. In 
this and other figures To is the condensation temperature (see 
text). 
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Fig. 3. The variation of To (in energy units of t=l) with A'^ 
for bosons in a Id lattice of size 1000 in a Id harmonic trap 
with k = 0.02 (top), 0.01 (middle), and for free bosons in a 
harmonic trap for k = 0.01 (bottom). 
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Fig. 4. Temperature dependence of the specific heat per par- 
ticle for 600 bosons. The dashed line is for free bosons in a 
harmonic potential. Solid lines are for bosons in a Id optical 
lattice of size 1000 in a harmonic potential for (from top to 
bottom): k = 0.001 (top), 0.005, 0.01, and 0.02 (bottom). 
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Fig. 5. Size effect in 2d for N = 540. The dashed curves are 
for k = 0.01 and lattice sizes: 30 x 30 (top), 40 x 40, 50 x 50, 
60 X 60, and 80 x 80 (bottom). Solid curves are for lattice 
bosons without harmonic potential for lattice sizes: 30 x 30 
(top), 40 x 40, 50 X 50, 60 x 60, and 80 x 80 (bottom). ksT is 
measured in energy units of t = 1. 
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Fig. 6. The variation of To (in energy units of t =1) with 
number of bosons. Solid lines are for an 80 x 80 lattice and 
with k = 0.05 (top), 0.01 (middle), and fc = (bottom). The 
dots are for bosons with only the harmonic trap potential with 
k = 0.01. 
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Fig. 7. Lattice size effect on specific heat in 2d for 3840 bosons Fig. 9. Specific heat per boson in 2d for 540 bosons in a 

in a harmonic trap with k = 0.01. The curves are for lattice 150 x 150 lattice in a harmonic potential with k — 0.01 (solid 

sizes: 150 x 150 (sohd fine), 100 x 100 (long dashes fine), 80 fine), k = 0.02 (dashed fine), and k = 0.05 (dash-dot line). 

X 80 (short dashed line), 50 x 50 (dotted line), and 30 x 30 The dotted curve (top) is for bosons in 2d harmonic trap and 



(dash-dot line). 



is independent of k. 
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Fig. 8. Specific heat per boson in 2d for 3840 bosons in a Fig. 10. NormaUzed DOS per site versus energy for an 100 x 
100 x 100 lattice in a harmonic potential with k = 0.01 (solid 100 lattice for k =0.01 (solid line) , 0.02 (dash-dot fine), and 



line), k = 0.02 (dash line), and k = 0.05 (dash-dot line). 



0.03 (dashed Une). 
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Fig. 11. Temperature dependence of condensate fraction for 

2500 bosons in a cubic lattice with 50 x 50 x 50 lattice sites and 
harmonic potential with k = 0.05 (top), k = 0.02 (dashed), k 
= 0.01 (dash-dot), and = (bottom). 
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Fig. 12. Bosons number dependence of To for bosons in a cubic 
50 X 50 X 50 lattice with k = 0.05 (top), k = 0.02 (dashed), k 
= 0.01 (dash-dot), and A: = (bottom). 
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Fig. 13. The effect of harmonic potential strength on the spe- 
cific heat per boson of 2500 bosons in a 50 x 50 x 50 cubic 
lattice. The results shown axe for k = 0.01 (solid line), k = 
0.02 (dashed line), and k = 0.05 (dash-dot line). Dotted line 
with a sharp peak is for 2500 bosons in a 3d harmonic trap. 
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Fig. 14. Specific heat per fermion in a Id optical lattice (of size 
3000 sites) in a harmonic potential with k = 0.0001 (dashed 
Une), 0.0005 (dash-dot), 0.01 (solid line). The filled circles are 
results in the absence of harmonic potential (i.e., for A; = 0). 
The number of fermions (iV) considered is 150. Tf is the Fermi 
temperature. 
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Fig. 15. Specific heat per fermion in a 2d optical lattice (of size 
250 X 250) in a harmonic potential with k = 0.01 (solid line), 
and 0.02 (dashed line). The dots are results in the absence of 
harmonic potential (i.e., for A; = 0). The number of fermions 
considered is 600. 




Fig. 16. Specific heat per fermion in a 3d optical lattice (of 
size 100 X 100 x 100) in a harmonic potential with k = 0.01 
(solid line), 0.02 (dashed line), and 0.05 (dash-dot line). The 
dots are results in the absence of harmonic potential (i.e., for 
k = 0). The number of fermions considered is 1000. 



